Midinfrared Third Harmonic Generation from Macroscopically Aligned 
Ultralong Single- Wall Carbon Nanotubes 
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We report the observation of strong tliird liarmonic generation from a macroscopic array of aligned 
ultralong single-wall carbon nanotubes (SWCNTs) with intense midinfrared radiation. Through 
power-dependent experiments, we determined the absolute value of the third-order nonlinear optical 
susceptibility, x'^' , of our SWCNT film to be 5.53 x 10"!'^ esu, three orders of magnitude larger than 
that of the fused silica reference we used. Taking account of the filling factor of 8.75% for our SWCNT 
film, we estimate a x'^' of 6.32 x 10~^^ esu for a fully dense film. Furthermore, through polarization- 
determining the 
component. 



dependent experiments, we extracted all the nonzero elements of the x tensor 
magnitude of the weaker tensor elements to be ~l/6 of that of the dominant x 
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PACS numbers: 42.65.-k, 42.65.Ky, 42.70.Km, 78.67.Ch 

Carbon nanomaterials, i.e., carbon nanotubes and 
graphene, attract much attention both from fundamen- 
tal and apphed viewpoints. These novel low-dimensional 
systems possess unique band structure and extraordi- 
nary properties that are promising for a variety of 
applications^^ Single-wall carbon nanotubes (SWC- 
NTs), in particular, are ideal one-dimensional systems 
for basic optical studies as well as for multi-wavelength 
photonic devices due to their diameter-dependent, direct 
band gaps. - Optical properties of SWCNTs have been ex- 
tensively studied during the last decade, and much basic 
knowledge has been accumulated on how light emission, 
scattering, and absorption occur in the realm of linear 
opticsiii^i^ However, nonlinear optical properties of car- 
bon nanomaterials remain largely unexplored although a 
number of interesting predictions existi^^— 

Theoretical calculations predict large non-resonant 
third-order nonlinear optical susceptibilities, ^^'^ ~ 10~^ 
to 10^^ esu, for SWCNTs, which vary rapidly with the 
tube diameterJ^"— Measurements of non-resonant y^"^^ 
of SWCNTs have been performed using some third- 
order nonlinear optical processes, including four-wave 
mixingiSii 



l^il^ and nonlinear refraction and absorption^"— 
where x^^^ was measured to be on the order of 10"^° to 
10"!^ esu. De Dominicis et al^ observed third harmonic 
generation (THG) from SWCNT films using nanosecond 
pulses of 1064 nm radiation, but the absolute value of 
X^'^-' was not quoted. 

Here, we have made the determination of the value of 
the x^^^ for SWCNTs via THG, by comparing our third 
harmonic intensity to the intensity of a reference material 
with well-known x^'^''- We have also made the determina- 
tion of the strengths of all the nonzero tensor components 
of the x*-'^^ tensor, by performing polarization dependence 
experiments on our highly aligned SWCNT sample. 



We measured THG from highly-aligned SWCNT films 
were fabricated via the process described in Ref. i2Sf To 
make the sample, carbon nanotube arrays were grown 
vertically via chemical vapor deposition on a Fe catalyst- 
lined substrate. The lines of catalysts were separated by 
a distance of 50 /im, and the self-supporting carbon nan- 
otube arrays were grown to a height specified by growth 
time. The vertically aligned carpet was then separated 
from the catalyst substrate via post-growth H2/H2O va- 
por etch to release chemical bonds between catalyst par- 
ticles and the nanotubes and then deposited horizontally 
onto a sapphire substrate. This process resulted in a 
large-area thin film of long, extremely well-aligned car- 
bon nanotubes on sapphire, with a small amount of over- 
lap. The film thickness was measured via vertical scan- 
ning interferometry to be 1.6 /xm. The films acted as 
nearly perfect polarizers for terahertz light^"— indicat- 
ing the high degree of alignment of the SWCNTs. 

THG measurements were performed using linearly- 
polarized midinfrared (MIR) radiation from an optical 
parametric amplifier pumped by a Ti:Sapphire-based 
chirped pulse amplifier (CPA-2010, Clark-MXR, Inc.). 
The MIR pulses generated had a wavelength of 2.1 /im, 
a repetition rate of 1 kHz, a pulse- width of ~300 fs, and 
a pulse energy of ^10/iJ. The MIR beam was incident 
on a CaF2 window that acted as a 97:3 beam splitter, 
where the reflected light was used to measure the funda- 
mental intensity with a liquid nitrogen cooled mercury 
cadmium telluride detector. The transmitted radiation 
was focused down to a spot size of ^100 /im in diam- 
eter, which allowed us to attain pump fluences of up 
to lOOmJ/cm^. Our sample was placed at normal inci- 
dence at the focal position, and a third harmonic signal 
at 700 nm was produced. The third harmonic beam was 
filtered by a monochromator and measured by a photo- 
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FIG. 1: (color online) (a) Third harmonic spectrum gener- 
ated at 700 nm (left, black) from highly aligned SWCNTs 
on sapphire, and the fundamental spectrum (2100 nm, right, 
red), (b) Shift in third harmonic wavelength (left, black) from 
700 to 500 nm due to a shift in the fundamental from 2.1 to 
1.5 /xm (right, red). The SWCNTs are ahgned parallel to the 
incident fundamental, and the induced third harmonic is po- 
larized parallel to the fundamental. 
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multiplier tube. In power-dependent experiments, the 
fundamental passed through a variable neutral density 
filter, to compare the measured third harmonic signal 
versus the incident fundamental power. In polarization- 
dependent experiments, the third harmonic signal was 
passed through a rotatable polarizer such that the mea- 
sured signal was either perpendicular or parallel to the 
incident fundamental polarization. The sample was then 
rotated through an angle (f) about its normal, to deter- 
mine the nonzero tensor components that contribute to 
the overall third harmonic signal. 

Figure [Ija) shows a spectrally resolved fundamental 
signal at 2.1 fim and its third harmonic signal at 700 nm 
generated from our SWCNT sample. Changing the fun- 
damental input wavelength caused a subsequent shift in 
the third harmonic signal, as shown in Fig. [ijb). The 
bulk sapphire substrate did not produce a measurable 
third harmonic signal, and thus, the total measured sig- 
nal can be attributed to the nanotube film. 

Figure [5] shows that the intensity of the third harmonic 
produced by the SWCNT sample and fused silica refer- 



FIG. 2: (color online) Power dependence of the 700 nm third 
harmonic versus the 2100 nm fundamental for the (a) 1.6 fim- 
thick highly-aligned SWCNT sample and the (b) 26.5 /xm- 
thick fused silica reference. Third harmonic shows cubic 
power dependence with the incident fundamental, as shown 
by the slope of the log/log plot (red). 



ence varies with the cube of the fundamental, as shown 
by the straight line power fits on this log/log plot. It can 
be seen that the same intensity produced by the carbon 
nanotube sample is achieved at an order of magnitude 
less input fundamental power. This order of magnitude 
difference between the fused silica and the carbon nan- 
otube cases shows that the x^^^ for the SWCNT sample is 
about two orders of magnitude larger than that for fused 
silica, as suggested by the coefficients of the power fits. 

In order to extract x*^^^ values for both samples from 
our experimental data, we need to consider phase match- 
ing conditions. The intensity of the third harmonic gen- 
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erated in a film of thickness L is given by 
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where /(^ (A^j) is the intensity (wavelength) of the funda- 
mental, Ak = k^i^ — Ski^ is the phase mismatch between 
the third harmonic and the fundamental, and Uui and 
nsui refer to the index of refraction at the fundamental 
and third harmonic frequencies, respectively. The so- 
called coherence length Lc = 27r/|Afc| = |n3^ — 
is a measure of the distance over which the fundamental 
and third harmonic remain in phase. Because x^'^^ is a 
fourth-rank tensor, the i*'' polarization component of the 
generated third harmonic field after traveling through a 
medium of length L is defined as 
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where the quantity in the bracket is the nonlinear di- 
electric polarization of the third harmonic, Ei^j, E^^ki 
and Ei^_i are the fundamental field in the j, /c, and I di- 
rections, and AcNT is a constant containing information 
about the fundamental-third harmonic interaction: 
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Because the SWCNT film is made of mostly air (as 
shown below), we assume the index of refraction to be 
nearly 1 for both the fundamental and third harmonic fre- 
quencies. Further, the sample thickness L is quite small 
compared to the expected coherence length. Thus, we 
take the the phase mismatch for the third harmonic and 
fundamental to be negligible: AkL <C tt/2. 

The thickness of the fused silica reference is 26.5 /im 
(measured via Fabry-Perot fringes using Fourier trans- 
form infrared spectroscopy), which is within the cal- 
culated group- velocity walk-off length assuming 300 fs 
pulses and on the order of the coherence length, Lc = 
17 fim. Because of the finite phase mismatch, the third 
harmonic intensity produced by the fused silica refer- 
ence sample depends somewhat sensitively on the sample 
thickness. The measured intensity for our fused silica ref- 
erence is less than the signal that would be produced if 
the sample thickness were exactly equal to the coherence 
length of 17 /im. This must be taken into account since 
the quoted value of x^'^'' for fused silica was measured at 
a multiple of the coherence length<^ 

By taking the ratio of the third-harmonic intensities 
of the SWCNT sample and the fused silica reference and 
solving for Xcnt' obtain 
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Because each of the power-dependence plots [Fig. [U^a) 
and Fig. [2Kb)] can be fit with a power function. 
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we compare the coefficients, C, to calculate the value of 
x'-'^' for the SWCNT film, using the relationship 
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By inserting all of our known quantities into Eq. Q, we 
find that the absolute value of x*^^ for the SWCNT film 
is 5.53 X 10~^^ esu. However, it is important to note that 
the SWCNT film is not fully dense. The carbon density 
of our film was measured^'^ to be 60mg/cm'^, with which 
we can calculate the filling factor to be 8.75%. Taking 
into account this filling factor, we find the value of x''^'' 
of a fully dense film to be 6.32 x 10^^^ esu, which is 
extremely high as non-resonant x'''^'' for any material. 

Far from resonances, one can make an order-of- 
magnitude estimate for x'-'^'' as 
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where /i — ej/{Ec — Ey) ~ evpjuj is the dipole matrix 
element of the interband optical transition'^* between va- 
lence and conduction band states of energies E^ and i?c, 
7 = (^3/2)070 = hvF, a = 2.46 A, 70 = 2.89 eV is 
the transfer integral, vp — c/300 is the Fermi velocity of 
graphene. An order-of-magnitude estimate of x'-"'^-' in the 
effective-mass kP description'^^ can be obtained by 
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where summation is performed over all one-dimensional 
electron fc-states, A = ttR^, Rt is the nanotube radius, 
g = A is the total dege neracy of an electron /c-state, 
Ec.vik) = 7v/K(^iF+fc2, Kin) = (l/i?t)(n- v/3), and 
the quantum numbers n and i' are defined in Ref. Issl . 
When far from any resonances and van Hove singular- 
ities at fc = 0, one can replace the summation by a 
typical wave number involved in the optical transition, 
fc/(27r) ~ hu! / (ATT-f) , and take ^ evp/uj. Assuming the 
difference in the occupation numbers /„ — fc to be equal 
to 1, taking the average radius R = 3 nm,— and using 
the frequency corresponding to the 2.1 /im wavelength, 
we arrive at x^^^ ~ M^/ (t^A'^) and x^'^^ ~ 9.8 x 10~^^ esu, 
in reasonable agreement with the measured value. 

The magnitude of x'''^' for SWCNTs can be also es- 
timated from what is known for graphene. The x'^'* 
of graphene at near-infrared wavelengths has been mea- 
sured to be ~10~'' esur^ which is the "bulk" suscepti- 
bility of a "graphene material" obtained by dividing the 
measured two-dimensional (sheet) susceptibility by the 
thickness, dg, of monolayer graphene. Using Eq. ^ with 
T^(i) ^ e^/ihwdg) and the value of w corresponding to 
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FIG. 3: (color online) Experimentally measured and theoret- 
ically calculated angular dependence for a THG signal po- 
larized (a) parallel and (b) perpendicular to the fundamen- 
tal, considering the x''^' tensor contribution relationship is 
aXzzzz ~ Xzzxx- The theoretical fits (red) show depen- 
dence for aRil/6. 4> = Q corresponds to light polarization 
parallel to the nanotube axis. 



our pump wavelength (2.1 /im), we obtain x^^^ ~ 3 x 
10~® esu, similar to the bulk susceptibility mentioned 
above. The effective x'^^ of a "SWCNT material" can 
then be obtained by "rolling up" graphene, i.e., diluting 
the volume by R^/dg ^ 150, which leads to the same 
ballpark estimate of x*-"^^ ^ 10^^° esu. 

Since our SWCNT film contains macroscopically 
aligned ultralong nanotubes, we can also determine the 
components of the x^^^ tensor. The measured third har- 
monic signal had both parallel and perpendicular polar- 
ization components relative to the incident fundamental 
polarization, as shown in Figs. [Sja) and ISj^b), respec- 
tively. There is strong (j) dependence in the third har- 
monic signal both when the third harmonic is polarized 
parallel and perpendicular to the fundamental. When the 
third harmonic is polarized parallel to the fundamental 
[Fig. [IJa)], the intensity measured when the nanotube 
axis is parallel to the fundamental polarization (0 — 0° 
and 180°) is almost two orders of magnitude larger than 



that of the case where the nanotube axis is perpendicu- 
lar to the fundamental (0 = 90° and 270°). When the 
measured third harmonic is polarized perpendicular to 
the fundamental, there are four peaks and valleys as the 
nanotube direction (f) is scanned. In both the parallel 
and perpendicular cases, the measured third harmonic is 
nearly zero when the fundamental polarization is perpen- 
dicular to the nanotube axis. 

These polarization-dependent results give us signifi- 
cant insight into the relevant nonzero x^'^'' tensor ele- 
ments for SWCNTs. Here we use Eq. 1^ to analyze 
these data. The incident fundamental is polarized paral- 
lel to the nanotube axis (z axis) when (f> = 0°. When the 
sample is rotated through an angle (p about its normal, 
we project this fundamental onto the axial (z) and radial 
(a;) directions of the carbon nanotubes as 

^ E^icz cos (j)- e.j, sin <j>). (9) 

By symmetry, the nonzero components of the x^^^ ten- 
sor for a SWCNT are required to obey the following 
relationships22i2^ 

0,^(3) ^ (3) ^ (3) ^ (3) ^ (3) ^ (3) ^ (3) 

^A.zzzz A,zzyy A.zyzy A.zyyz A,zzxx A^zxzx A^zxx 

(10) 

where Xzzzz = Xcnt ^^'^ a (0 < a < 1) is the ratio 
of the weaker tensor components to the dominant tensor 

(3) 

component, Xzzzz- Thus, the third harmonic field along 
the nanotube axis is 

Eauj^z = ^cNTXcNT^^(cos^<^ + 3acos0sin2 0). (11) 

After projecting the nanotube coordinate system back to 
the coordinate system of the incident fundamental, we 
find that the component of the induced third harmonic 
field parallel to the incident fundamental is 

E3uj,\\ = v4cNTXcNT-^<iJ(cos^'?^ + 3acos(/)sin^(/))cos0, 

(12) 

and the corresponding intensity is 

hu^.w = |^cNTnXcNTl'^'(cos^'^ + 3acos2 0sin2 0)2. 

(13) 

Simulations show that the peak that would arise purely 

f3) 

from the dominant tensor component Xzzzz becomes 
broadened as the contribution from the weaker tensor 
components increases (i.e., as a increases from zero). 
Similarly, the induced third harmonic intensity polarized 
perpendicular to the incident fundamental is 

^31.,;, = |AcNTnXcNTp-^<2(cos^ '/'Sin0+3acos(/)sin^ (p)^. 

(14) 

As the contributions from the weaker tensor components 
increase (i.e., as a increases from zero), the peaks of the 
induced third harmonic shift and sharpen. 

We fit our data with the theoretically determined fit- 
ting functions in Eqs. (IT51) and ([Ti| . By allowing the 
value of a to be the only adjustable parameter, we found 
that the measured third harmonic signals (black line with 
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cross markers) are in excellent agreement with the theo- 
retically calculated fits (red line), as shown in Figs.[3Ka) 
and El^b) . Not only do the fits correlate extremely well 
with the measurements, but the fit parameter a is ap- 
proximately equal to 1/6 in both cases. This value indi- 
cates that the dominant x*-"^) tensor component, xSzz, is 
six times larger than the weaker components. 

In Fig. [3]^ a) , the third harmonic signal is small but fi- 
nite when the SWCNTs are oriented perpendicular to the 
incident fundamental {(p = 90° and 270°), whereas the- 
ory based on either intraband or interband transitions 
predicts zero third harmonic signal at these points. Con- 
tributions from cross-polarized transitions, e.g., £"„ to 
En±i, could be a source of this finite signal in the perpen- 
dicular configuration although such transitions are weak- 
ened by the depolarization effect 

In summary, we successfully observed third harmonic 



generation in highly aligned SWCNTs on a sapphire 
substrate. Through power-dependent measurements, we 
were able to determine x*''^'' to be 5.53 x 10^^^ esu for 
the film. With an estimate filling factor of 8.75% for our 
film, x^"^) value of a fully dense film would be 6.32 x 
10^^^ esu. Furthermore, through orientation-dependent 
experiments we were successfully able to extract all the 
relevant nonzero x'^'* tensor elements. We found that 
weaker x^'^'* tensor elements are approximately 1/6 the 

(3) 

strength of the dominant Xzzzz- 
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